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We show that a normal (single particle) current density J x transverse to the ordering wavevector 
2kpz of a charge density wave (CDW) has dramatic effects both above and below the CDW depinning 
transition. It exponentially (in J x ) enhances CDW correlations, and exponentially suppresses the 
longitudinal depinning field. The intermediate longitudinal I-V relation also changes, acquiring a 
linear regime. We propose a novel "current effect transistor" whose CDW channel is turned on by a 
transverse current. Our results also have important implications for the recently proposed "striped 
phase" of the high-T c superconductors. 



PACS: 7f.45.Lr, 64.60.Ht 



While the equilibrium properties of charge density 
waves (as well as other disordered periodic media) 
are by now fairly well understood, recent attention has 
focussed on their far less well understood non- equilibrium 
dynamical properties. 

All classical treatments of CDW's to date have in- 
cluded only the local displacement u(r, t) of the CDW, 
which is related to the electron density via p(r, t) — 
Pa + piCos(2fci?(z + u(r,t)), as the only important hy- 
drodynamic mode of the problem. 

In this Letter we will argue that while such an approach 
is correct for the statics, the total electron density p is 
crucial for the correct description of the non- equilibrium 
dynamics. We will present a dynamical model that in- 
corporates this additional mode and explore its conse- 
quences. 

This additional mode allows current to flow even when 
the CDW itself is stationary. The effects are particu- 
larly dramatic when the current flows perpendicular to 
qo = 2kpz, the ordering wavevector of the CDW. Once 
this transverse current J x exceeds a crossover value J c , 
it makes the CDW much more ordered. The correlation 
length of a three dimensional (3d) CDW along the direc- 
tion (x) of the transverse current obeys 



{Jx 



,2{J X -J C )/J C 



for J x < J c 
for J x > J c 



(1) 



where £l is the "Larkin length" (see below) of the CDW 
in zero transverse current, which is finite due to the ran- 
dom pinning of the CDW by impurities. 

For directions perpendicular to the transverse current 
the correlation lengths also grow exponentially with J x , 
but with precisely half the growth rate: 



The "critical" current J c in the above expressions is: 



Jc 



,E T {0){k F £ L )(p n /p cdw ) x 0(1) 



(3) 



where Et(J x = 0) is the threshold electric field along 
2k f z necessary to depin the CDW, and CToo the conduc- 
tivity in that direction for E z ^S> Et(0). For typical 
NbSe3 samples J c ~ 10 3 — 10 4 Amps/cm 2 . These quite 
high values can be greatly reduced by increasing disorder 
(lowering £ L ). 

In addition to being directly observable by X-ray scat- 
tering measurements of the CDW correlation length, 
these exponential dependences of the correlation lengths 
on J x also have striking consequences for the I-V char- 
acteristics along z. If we apply a "longitudinal" electric 
field E z along z while maintaining a fixed J x , we find 
the (zero temperature) threshold depinning field Et(J x ) 
required to first make the CDW move is given by 



Et(J x ) 



Et(0) , 
Et(0) 



3/2 



-2(J X -J c )/Je 



for J x < J c , 
for J x > J c . 

(4) 



for J T > J r 



(2) 



Correspondingly, the shape of the high velocity part 
of the I-V is also drastically modified. We find that 
for longitudinal fields E z in the range, Et{J x ) « 
E z « E T {0)(J x /J c ) 2 , the I-V is Ohmic (i.e. linear) 
with the longitudinal conductivity given by o~ z (J x ) = 
CToo (1 - (Jc/Jx) x 0(1)). For the largest fields E z » 
Et(0)(J x / Jc) 2 , Jz(E z ) crosses over to the J x — re- 
sult: J Z {E Z ) = a^E, (l - cy/Er{0)/E z ) . § All of our 
results for the longitudinal I-V characteristics in the pres- 
ence of a transverse current are summarized in Fig.l. 

Our results suggest a novel "current effect transistor" 
(CET) , in which the threshold field along qo is controlled 
by a current driven transversely to qo. 
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All of the above results should also apply to the re- 
cently proposed "striped phase" of the high-T c cuprate 
oxide superconductors, H since the macroscopic sym- 
metries of that system (i.e., unidirectional, lattice- 
incommensurate modulation of the charge density) are 
the same as those of a CDW, as is the ubiquitous pres- 
ence of disorder. 
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-2(Jx- Jc)/J, 



Fig.l:(a) Schematic I-V curve illustrating new intermedi- 
ate Ohmic high velocity regime that exists for transverse 
currents J x > J c . 

We now present our model and derive these results. 
General principles dictate that there are two types of long 
wavelength hydrodynamic variables: Goldstone modes 
associated with broken symmetries, and conserved vari- 
ables. In CDW's, the phonon mode u associated with the 
breaking of translational symmetry is the sole Goldstone 
mode, while the total electron density p is a conserved 
variable. Contrary to the naive expectation, p is not 
determined solely by the compression of the CDW (i.e. 
dp =/= —pcdwd z u). This is in direct analogy with ordinary 
3d crystals, for which the density of vacancies and inter- 
stitials (or equivalently the total particle density) must 
be included in the proper hydrodynamic description 

Although the total charge density p is a conserved field, 
in CDW's, because of long-range Coulomb interactions 
||, it does not constitute a slow hydrodynamic mode, 
unlike neutral systems, e.g. the transverse smectic phase 
of a driven vortex lattice || . 

In contrast, in driven CDW's one must consider to- 
tal charge density dynamics to even obtain the correct 
equation of motion, to which we now turn. 

We consider a pinned CDW with an ordering wavevec- 
tor qo = 2kp2 and explore the consequences of a "nor- 
mal" current J. Our main conceptual point is that, in 
contrast to all previous treatments such a current makes 
even the stationary CDW a highly non- equilibrium sys- 
tem in which the J — ► —J symmetry is broken. Lacking 



the spatial inversion symmetry along the direction of such 
a single-particle current, general symmetry principles 
therefore dictate that the equation of motion for u(r,t) 
must admit and therefore will contain non-equilibrium 
terms that break this inversion symmetry. The most im- 
portant of these is the "convective" term v ■ Vm. The 
striking effects of this should be experimentally observ- 
able. 

We explore these in a transverse geometry, in which the 
charge current J flows transversely to the CDW ordering 
wavevector 2/cfz, i.e. J = J x 5c. The hydrodynamics of 
a CDW in the presence of such transverse current J x is 
described by 

j(dt + v ■ V)u(r, t) = KV 2 u + F(r, u) + ep cdw E z , (5) 

where E z is the externally imposed electric field along 
\iF and p c dw is the average CDW electron number den- 
sity. We take the quenched random force F(r,u) in 
Eq.|5| to be Ga ussian distributed with zero mean and 
F(r, u)F(r', u') — A(u — u')S(r — r'), where the overbar 
denotes a disorder average. The function A(u) is periodic 
with the CDW lattice spacing, a = 2ir/2kF- 

Because Galilean symmetry is broken by the underly- 
ing lattice and quenched disorder the "convective" cou- 
pling v is different from (but probably on the order of) 
the "normal" electron velocity v = J / p n e, where p n is 
the "normal" electron density. For simplicity of notation 
we have taken K x = K y = K z = K , however all of our 
results can be trivially extended to anisotropic elasticity. 

The essential difference between our model Eq.[| and 
the commonly used Fukuyama-Lee-Rice (FLR) model 
is the "convective" v x d x u term, which physically arises 
because a tipped (d x u ^ 0) "layer" of the CDW de- 
flects the transverse normal current downward by an an- 
gle 9 ~ d x u, leading to a reaction force back on the CDW, 
proportional to d x u and the normal current J x . 

This "convective" term leads to a velocity-dependent 
crossover length scale L c (v x ) — K/(^v x ). On length 
scales L x (along x) < L c , the convective term is unimpor- 
tant and the CDW behaves as though it is in equilibrium; 
|]l0f for L x > L c , non-equilibrium effects set in. 

In spatial dimensions d < 4, even arbitrarily weak pin- 
ning destroys the translational order of the CDW for 
length scales longer than the so-called Larkin length scale 
£l . (lljnj On this length scale the disordering effect of the 
random pinning on u just balances the "elastic" forces 
that try to keep the CDW ordered. 

Clearly, for vanishing or small transverse currents 
(small v x ) the standard FLR-Larkin (FLRL) calcula- 
tion of £l, which ignores v x , applies, giving in 3d £l = 
K 2 a 2 /A(0). 

Given the two length scales £l and L c (v x ), one can de- 
fine a crossover transverse velocity v c by £l = L c (v c ), 
giving v c = K/{^ L ) = A(0)/( 7 Xa 2 ). Taking v x sa 
Jx/{Pne), we can obtain a crossover current density 
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Jc = Pn^K/ Using the standard FLR result for the 
longitudinal threshold electric field Et — Kaj '(pcdwe^l) 
and the large-field limit of the longitudinal CDW conduc- 
tivity (Too = e2 Pcdw/li we obtain Eq|| with J c expressed 
in terms of experimental observables. For J x < J c , 
the equilibrium FLRL calculation of £x, is valid and the 
CDW translational correlation length is independent of 
the transverse current J x . 

Interesting transverse current dependent phenomena 
occur for J x > J c . In this regime, we can calcu- 
late the new CDW correlation length £,±(v x ) in the _L- 
direction by asking how big in the direction perpen- 
dicular to x the system would have to be before the 
mean-squared fluctuations (u 2 ) in the position of the 



CDW exceed 



For a finite system of width L±, in 



d < 4 dimensions, the mean-square phonon fluctuations 
(w 2 (r)) = J q±>L -i dq x d d ~ 1 q±{\u(q)\ 2 ) can be calculated 
by spatially Fourier transforming Eq.^ and looking for 
the static solution for the spatial Fourier transform field 
u(q). 

For a system with Lj_ < £±(v x ), the u fluctuations re- 
main smaller than the lattice spacing a, and so we can 
replace F(r, u) with its u — value. The static equation 
of motion then reads 



(ijv x q x + Kq 2 )u(q) = F(q; u = 0) , 



(6) 



where -F(q; u = 0) is the spatial Fourier transform of 
F(r;u = 0). Solving this for u(q), and then computing 
(|u(q)| 2 ) from the result, we obtain 

A(0) 



(Kq)l : 



(7) 



Inserting this into expression for (u(r) 2 } and defining 
£,±(v x ) as the value of L± at which (u(r) 2 ) = a 2 we obtain 
an implicit equation for £±(v x ) 



4-rf 



(8) 



where f(y) = 2(4-d) dxx^/ig+g^+g^g^), with 
9±(x,y) = x 2 +y 2 /2±y{x 2 + y 2 /4) 1 / 2 . 

Evaluating the integral in d = 3 in the J x << J c and 
J x >> J c limits and using v x sa Jx/(Pn e ) gives Eq.^ 
for £,±(J X ), which matches the J x = Larkin length at 
J x = J c , i.e., £,±(J C ) « £l- The correlation length £ X (J X ) 
along the transverse current J x (v x ) is easily computed 
from £_l(J x ) to be given by Eq.0. 

Physically, the exponential dependence of the CDW 
correlation length for J x > J c arises due the strong 
suppression of the CDW roughness by the transversely 
moving normal charge carriers, which "stiffen" the CDW 
elasticity by momentum transfer with it. The transverse 
convective term lowers the upper critical dimension (be- 
low which the random pinning disorders the CDW) from 
d uc {v x < v c ) — 4 to d uc (v x > v c ) = 3. 



The strong exponential dependence of the Larkin 
length on a transverse single particle current J x has a 
number of striking consequences which should be exper- 
imentally observable (aside from being directly probed 
in X-ray scattering). It is easy to show that at low 
temperatures and weak disorder, the threshold longitu- 
dinal field (along the ordering wavevector 2kpz) is given 
by E T (J X ) = E T (0) (V l (0)/V l (J x )) 1/2 , where E T (0) = 
Ka/(ep c dw£,j j ) is the J x = longitudinal threshold elec- 
tric field and Vl(J x ) is the transverse current dependent 
Larkin volume. This decay of the longitudinal thresh- 
old field Et with Larkin volume is physically associated 
with the V^V statistic reduction in the strength of the 
random pinning force, when averaged over the correla- 
tion volume V(J X ) = £,j_(J x )£x(J x )- Using our results for 
£ x ,x(J x )i Eq.[yfJ in the expression for Et(J x ) above, we 
obtain Eq.||. 

We now turn to the effects of the transverse current 
J. x for E z well above the longitudinal threshold field 
Et{Jx)- Perturbing about the uniformly moving CDW 
state uq = v z t, we seek a solution to EqJ^ of the form 
u (r, t) = v z t + Su (r, t), where v z is the mean velocity of 
the CDW, to be determined self-consistently, and Su (r, t) 
a zero mean fluctuation (assumed small) about this uni- 
formly moving state. Inserting this Ansatz into Eq||, 
expanding for small Su, and averaging over the disorder, 
we have v z — r )~ 1 ep c d w E z + Sv z , with 



-ySv z rj d u F(r,u)Su(r,t) 



(9) 



where we have used the facts that F(r, v z t) — ~S/ 2 5u = 0. 
We can obtain d u F8u from the equation of motion for 
Su(r,t). To leading order in 5u, that equation reads: 



(Tfdt + -fo ■ V - K V 2 )8u = Y, F(r, q n ) e lq " v * 



(10) 



where we have used the fact that F(r, u) is a periodic 
function of u and q n = nqo, with n € Z. 

Solving Eq.[l0| for Su and using it inside Eq.^ for Sv z 
we obtain 



lb~v z 



E 



d d q \ n q n (v z q n 



v ■ q) 



(27r)d 1 2 (v z q n ~v-q) 2 + (Kq 2 ) 2 



(11) 



We can verify a posteriori that the wavevector integral in 
this expression is dominated by q z -C q n . This together 
with v ~ v implies v ■ q ss v x q x , so: 



lSv z 



d d q 2A Q v x q (q 1 - q x ) 
{2nY 1 2 v 2 x (q 1 -q x ) 2 + K 2 q i 



(12) 



where we have defined q\ = qov z /v x and kept only n = 1; 
higher random force harmonics (n > 1) do not qualita- 
tively change the results. 

We can asymptotically evaluate the integral in Eq.|T^ 
in the limits: (i) q* <C qi and (ii) q* 3> qi, where 
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q* is the characteristic value of q x that dominates the 
integral. A simple analysis shows that in (i) q* = 
\Jl v z<lo/K. The condition q* < q ± then leads to the 
constraint E z /E T (0) > Max{l, (J x / J c ) 2 }. Ignoring q x 
relative to qi in Eq. 12 we find that the SCF result 



M holds at these highest longitudinal fields: Sv z /v z — 
-{E T {Q)/E z )^- d ^ 2 x 0(1). For J x < J c , this is the 
only regime. 

For J x > J c , we find that a new regime (regime (ii)) ex- 
ists, when the longitudinal fields E z /Et{0) < (J x /J c ) 2 - 
Here q\ is subdominant to q x and the integral above can 
be evaluated by a Taylor expansion in q± . Changing vari- 
ables q x — qi — > q x and expanding in qi we find 

d d q 



^ « ~SA q 2 Q K 2 [ p 
v z J (2?r) 



2 2 

TxT 



Y tfvlql + K 2 q±] 

The integral on the right hand side is finite in the infra- 
red and the ultraviolet for 1 < d < 4. All higher terms 
in this qi expansion are likewise finite. Evaluating the 
integral in d = 3 we obtain 

^-M x0(1) = _!i x0(1) 

v z K~/v x v x 

A more detailed analysis shows that the above 
expansion breaks down as Et{v x ) is approached from 
above. This can be seen simply by assuming that v z 
in Eq.|| is of order the v z we found above, and us- 
ing that result in a calculation of (Su(r,t) 2 ). We find 
that this approaches a 2 , invalidating our small 5u ex- 
pansion, for E z ~ the threshold field Et(v x ) found ear- 
lier. This breakdown signals the pinning of the CDW for 
E z < E T (v x ). 

This conclusion is consistent with our earlier observa- 
tion that J c is a crossover (as opposed to critical) trans- 
verse current, and with similar results for the related 
transverse smectic phase of a driven vortex lattice. || 
Our results for the transversely driven CDW I-V charac- 
teristics are summarized in Fig.l. 

Another observable quantity sensitive to transverse 
current is the Fourier transform I(r) of the structure 
function: 



over to the familiar super-universal (for 2 < d < 4) log- 
arithmic dependence on r [^2|. At scales larger than 
L c {v x ) the correlation function is also logarithmic (but 
with non- universal coefficient), now as a consequence of 
the fact that non-equilibrium effects have lowered the 
upper-critical dimension to d uc = 3. S For J x > J c 
the intermediate regime is absent and there is a direct 
crossover at the Larkin scale to the logarithmic growth 
associated with non-equilibrium effects. 

Finally we note that our results, in particular EqJJ, 
might have an interesting device application: the "cur- 
rent effect transistor" (CET). Imagine a channel made 
of e.g., NbSe3. |l3| At low T, and longitudinal fields be- 
low the threshold field Et, such a channel is insulat- 
ing. For E z > Et the CDW depins and the channel 
becomes highly conducting. Our analysis demonstrates 
that this longitudinal threshold field Et{J x ) is a strong 
function of the transverse current J x for J x > J Cl as 
summarized by Eq.|]. One can therefore turn on the 
CDW transport by a transverse current J x which depins 
the CDW by depressing the longitudinal threshold field 
(14) E T (J X ) ~ e -2(Jx-Jc)/Jc_ Though theoretically interest- 



(13) 



for J x < J c , 



where n( J x ) = 0(1) x J c /J x and gr)(z) is a scaling func- 
tion readily calculated to lowest order in weak disorder. 
It is important to keep in mind that the expression above 
is only valid for length scales smaller than the Larkin 
length. 

On scales larger than the Larkin length, the small u 
expansion of the random force is invalid. However, using 
the functional renormalization group developed in Ref. 
p|, we can show that u fluctuations are logarithmic for 
both J x < J c , and for J x > J c . In the former case, the 
phonon correlation function C (r_i_ , x) (which grows lin- 
early with r for r < £l ) first (at the Larkin scale) crosses 



ing, this "device" would almost certainly be too slow and 
low gain to be practical. 

In conclusion, we have argued that previous treatments 
of sub-threshold CDW dynamics have ignored an im- 
portant "permeation" (total charge density) mode. We 
demonstrate that this mode radically changes the CDW 
dynamics in the transverse current geometry. 
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